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The relic abundance and the scalar cross-section off nucleon for light neutralinos (of mass below 
about 45 GeV) are evaluated in an effective MSSM model without GUT-inspired relations among 
gaugino masses. It is shown that these neutralinos may provide a sizeable contribution to the matter 
density in the Universe and produce measurable effects in WIMP direct detection experiments. 
These properties are elucidated in terms of simple analytical arguments. 
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I. INTRODUCTION 

Most works on relic neutralinos consider supcrsymmetric schemes with a unification assumption for the gaugino 
masses Ms (i = 1, 2, 3) at the GUT scale Mqut ~ 10 16 GeV. This hypothesis implies that at lower scales the following 
relations hold: 

Mi : M 2 : M 3 = ai : a 2 : a 3 , (1) 

where the on [i = 1, 2, 3) are the coupling constants of the three Standard Model gauge groups. In particular, at the 
electroweak scale, Mew ~ 100 GeV, Mi and M 2 are related by the expression: 

Mi = - tan 2 9 W M 2 ^ 0.5 M 2 . (2) 

o 

However, there are theoretical arguments for considering supersymmetric schemes where the unification assumption 
on gaugino masses is not satisfied 

In the present paper we analyse properties of relic neutralinos in an effective Minimal Supcrsymmetric extension of 
the Standard Model (MSSM) where the GUT relation of Eq. (0) is relaxed. Previous papers where supersymmetric 
schemes without gaugino masses unification have been considered in connection with relic neutralinos include the ones 



;augmo masses unification nave been c 
reported in Refs. 0, S 01 El 0- B 0- EH EH E3 E3- Here we evaluate the neutralino relic abundance £l x h l and 
the neutralino-nucleon scalar cross-section o~^™^ , which is relevant to dark matter direct detection. In Sect ITT1 we 
define the supersymmetric scheme adopted in the present paper and in Sect. IHII we provide analytical considerations 
and numerical evaluations. Our conclusions are drawn in Sect. IIVI 



II. EFFECTIVE MSSM WITHOUT GAUGINO UNIFICATION 



We employ an effective MSSM scheme (cffMSSM) at the electroweak scale, defined in terms of a minimal number 
of parameters, only those necessary to shape the essentials of the theoretical structure of MSSM and of its particle 
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content. The assumptions that we impose at the electroweak scale are: a) all squark soft-mass parameters are taken 
degenerate: = nig] b) all slcpton soft-mass parameters are taken degenerate: mr = mr, c) all trilinear parameters 
are set to zero except those of the third family, which are defined in terms of a common dimensionless parameter A: 
A~ h = Ai = Arrig and Af = Amj. As a consequence, the supersymmetric parameter space consists of the following 
independent parameters: M 2 , H, tan/3, toa, f/if, mr, A and R = M1/M2. In the previous list of parameters we have 
denoted by fi the Higgs mixing mass parameter, by tan/3 the ratio of the two Higgs v.e.v.'s and by tua the mass of 
the CP-odd neutral Higgs boson. 

This scheme differs from the effMSSM which we employed for instance in Rcf. 01 m the f ac t thai we are relaxing 
here the gaugino unification relation, which was instead assumed in our previous works. The presence of the extra R 
parameter accounts for this fact. 

The neutralino is defined as the lowest-mass linear superposition of bino B, wino and of the two higgsino 

states HI, H°: 

X = a 1 B + a 2 W {3) +a 3 H° + ai H°. (3) 

Due to well-known properties of the neutralino and chargino mass matrices, one has that: a) for \i 3> M\,M 2 the 
neutralino mass is determined by the lightest gaugino mass parameter: m x ~ min(Mi, M 2 ), while the lightest chargino 
mass is set by M 2 : m x ± ~ M 2 (Mi does not enter the chargino mass matrix at tree-level); b) for /x <C Mi, M 2 both 
the neutralino and the chargino masses are primarily set by the Higgs mixing parameter: m x ~ \i ~ m x ± • 

LEP data put a stringent lower bound on the chargino mass: m x ± > 103 GcV, which converts into lower bounds 
on M 2 and fi: M 2 , fi > 103 GcV. This implies a lower bound on the neutralino mass of the order of about 50 GeV in 
the standard effMSSM, where the GUT relation of Eq.J2J) holds. On the contrary, the neutralino mass may be smaller 
when Mi <C M 2 , thus for small values of the parameter R. 

In the present paper we are interested in the phenomenology of light ncutralinos, therefore we consider values of 
R lower than its GUT value: Rgut — 0.5. For definiteness we will consider the range: 0.01 - 0.5. The ensuing light 
neutralinos have a dominant bino component; a deviation from a pure bino composition is mainly due to a mixture 
of B with HI, as will be shown in Scct lHI 51 

In our numerical analysis we have varied the MSSM parameters within the following ranges: 1 < tan/3 < 50, 
100 GeV < \n\, M 2 ,mg,mT < 1000 GeV, sign(» = -1,1, 90 GeV < m A < 1000 GeV, -3 < A < 3, for a sample of 
representative values of R in the range 0.01 < R < 0.5. This range for R, implemented with the experimental lower 
limit on M 2 of about 100 GcV, implies that the lower bound on the neutralino mass can be moved down to few GeV's 
for R ~ 0.01. 

We then implemented the following experimental constraints: accelerators data on supersymmetric and Higgs 
boson searches (CERN e + e~ collider LEP2 0] and Collider Detector CDF at Fermilab H3); measurements of the 
b — > s + 7 decay We wish to comment that the accelerator limits on the Higgs sector are taken into account 

by implementing the limits on the Higgs production cross-sections: e + e~ — > hZ and e + e~ — > hA (h and A are 
the lightest scalar and the pseudoscalar neutral Higgs bosons, respectively), which in turn imply a constraint on the 
coupling constants sin 2 (a — /3) and cos 2 (a — /3). Once these limits are applied, the absolute lower limit on the Higgs 
masses is m,A,mh ~ 90 GeV. The allowed light-Higgs mass range between 90 and 114 GeV is very often overlooked 
in studies of neutralino dark matter, where a flat limit of 114 GeV is applied to n ih- The light-Higgs mass range, 
even though difficult (but not impossible) to be achieved in SUGRA models IToL Eo J, i s nevertheless quite natural 
in the effMSSM and usually provides large detection rates for neutralino dark matter [19j. 

As for the constraint due to the muon anomalous magnetic moment a M = (g^ — 2)/2 we have used the interval 
— 160 < Aa M • 10 11 < 680, where Aa M is the deviation of the current world average of the experimental determinations 
(dominated by the measurements of Ref. jil]]) from the theoretical evaluation within the Standard Model: Aa M = 
a expt _ a SM _ rpjjg ran p- e we use f or /\ a i s a 2a interval, obtained by using for the lowest-order hadronic vacuum 
polarization contribution an average between the results derived from the e + — e~~ data |22j, |23j and from hadronic r 
decays |'22i . The Aa M constraint and the b — > s + 7 bound set stringent limits for the light neutralino sector of our 
models. 

Once also the relic abundance bound Q, x h 2 < 0.3 is applied (see Sect lHI Bj) in addition to the other experimental 
constraints discussed above, a lower limit of about 6 GeV is obtained for the neutralino mass in the class of models 
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with non-universal gaugino masses considered in this paper 1 . 



III. NEUTRALINO RELIC ABUNDANCE AND NEUTRALINO-NUCLEON CROSS-SECTION 

A. Some analytical properties for small m x 



The neutralino configurations which provide the highest values of direct detection rates are the ones dominated by 
(h, H) Higgs-exchange processes, which in turn require a gaugino-higgsino mixing. For these configurations, also the 
relic abundance is regulated by a (A)Higgs-cxchangc diagram in the \ — \ annihilation cross-section. 

Thus, to get an insight into the properties to be expected for our light ncutralinos we limit ourselves to the following 
approximate expressions, derived under the assumptions of Higgs-dominance and light neutralinos (notice however 
that full exact expressions both for the relic abundance VL x h 2 and for the neutralino-nucleon scalar cross section 



(nucleoli) 
^scalar 



are employed in the numerical evaluations to be discussed in the next Section). Under these hypotheses, 
the neutralino relic abundance is dominated by the s-wave annihilation in a bb pair (unless m x is very close to the 
fo-quark mass rrib, in which case the cc and tt channel are dominant): 
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and the elastic scattering cross section is: 



(nucleon) ^£.,2 2 
°scalar — ~~ iW Z TO rod 
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(4) 



(5) 



In the previous equations we have used the following notations: (a ann u)i n t is the integral from present temperature 
up to the freeze-out temperature of the thermally averaged product of the annihilation cross-section times the relative 
velocity of a pair of neutralinos; e is a quantity which enters in the relationship between the down-type fermion 
running masses and the corresponding Yukawa couplings (see, for instance, Refs. 0, |2f| and references quoted 
therein); m re d is the neutralino-nucleon reduced mass. The quantities Fh,H an d Ih,h are defined as follows: 



F,, 
F H 

Ih,H 



(— di sva.0\y + a 2 cos #w)( a 3 sin a + CI4 cos a) , 
(— ai smdw + 0-2 cos 0w)( a 3 cos ot — 04 sin a) , 



J2^' H m q {N\qq\N) 



(6) 



The matrix elements < N\qq\N > are meant over the nucleonic state. The values adopted here for m q < N\qq\N > 
are the ones denoted by set 1 in Rcf. [24|. We remind that uncertainties in the values of m q < N\qq\N > can give 
rise to an increase of the neutralino-nucleon cross section of about a factor of a few [25| . 
The angle a rotates H[ ^ and into h and H, and the coefficients k^' H are given by 



u— type 



= cos a I sin (3 . 



fcrf_ty pc = — sin a/ cos [3 — e cos(ct — /?) tan [3 . 
= sin a I sin/?, 



u— type 



k 



d-ty-pc = cos a/ cos /3 — esin(a — 0) tan/3 , 



(7) 



for the up-type and down-type quarks, respectively. 



1 This is at variance with the results of Ref. |33l . where a lower limit on the neutralino mass of 12 GeV has been deduced. Notice that 
the authors of Ref. l33l consider only the limit of very large m^, which strongly suppresses processes which involve A-exchange, in 
particular the neutralino annihilation cross section. On the contrary, we are considering also the light Higgs sector, which is effective in 
reducing the value of the neutralino relic abundance and therefore in allowing lighter neutralinos (see Ea. llUI I. 
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In the discussion which follows we only wish to establish some correlations implied by the dependence of il x h 2 and 

°f ^scala^ on t ne Higgs masses and the neutralino mass. For this purpose we rewrite the two previous expressions 
as follows: 

Q h 2 „ c [(2m x ) 2 -mj] 2 ^ 

rrn 2 M — rn^ I nm 2 11/2 ' 



(nucleon) ^ ~_ , q s 
"scalar — 4 ' W 
" l h 

with obvious definitions for C and D. Here mt stands generically for the mass of the one of the two CP-even neutral 
Higgs bosons which provides the dominant contribution to ^f°aiar ■ 

We now consider the case of very light neutralinos, i.e. m x <C \raA- Therefore we may further simplify Eq. JHJ) as 

n x h 2 ~ c 2 m ^ 2/ /2 (io) 

The largest neutralino-nucleon scattering cross sections occur when both m,h and toa are close to their experimental 
lower bound (m/, ~ ~ 90 — 100 GeV) and tan /3 is relatively large, in which case also the couplings of Eqs. i|6l7|) 
between neutralinos and down-type quarks through /i-exchange are sizeable 0|- I n this case, from Eqs. © and Q10|l 
one derives the range of o'i^ai & t fixed value of m x (always in the regime m x <C Im^): 



'scalar " u vcliLlc ui '"X ^""d" 3 10 6""= '"X ^ 2 

- ^(nucleon) ^ ■£ 

^min 



_ <^ ^nuLitun; , — fill 

m 2 [l_ m 2/ m 2]l/2 (fi x /,2 )max ~ Scalar 



where rrih,min stands for the experimental lower bound on m/j. The lower limit to ^""lar displayed in Eq. I|ll(l 
provides a stringent lower bound on ^gcaiar° ^ or ver y light neutralinos. This feature will show up in the numerical 
evaluations presented in the next Section. The upper bound on ^scalar is instead determined by the lower limit on 
the Higgs mass to/,. 

By the arguments given above, it turns out that in the small mass regime (m x <C ^tra) the upper bound on the 
relic abundance fl x h 2 < 0.3 establishes a constraint between the otherwise independent parameters m x and toa (see 
Eq.JTIJ). 



B. Numerical results 

We turn now to the presentation of our numerical results. In Figs llafJTbl we give the scatter plots of the quantity 
£ a scalar m terms of the neutralino mass for different values of the parameter R. The quantity £ is defined as the 
ratio of the local (solar neighbourhood) neutralino matter density to the total local dark matter density: £ = p x j p\ oc . 

In Figs llatiTbl we plot the quantity ^cr^^r > rather than simply &gc&lar ■> m or der to include in our considerations 
also neutralino configurations of low relic abundance (i.e. cosmologically subdominant neutralinos). We recall that, 
from experimental measurements of the direct detection rates, only the product ^scalar may be extracted, and 
not directly ^scaiar 50 '^ • The quantity £ is derived here from the relic abundance by the usual rescaling recipe: £ = 
min(l, Q x h 2 / [flh 2 ] m i n ) , where the minimal value of relic abundance which defines a neutralino as a dominant dark 
matter component has been fixed at the value f0fe 2 ] m in = 0.05. Q x h 2 and ^"aiar are evaluated according to the 
procedures and formulae described in Refs. |24j,|26|. 

Figs llaf)lbl displays quite remarkable properties of the light relic neutralinos from the point of view of their de- 
tectability by WIMP direct measurements. These properties are easily understandable in terms of the analytic 
arguments presented in the previous Section. For instance, in each panel at a fixed value of R < 0.1, there is a 

characteristic funnel pointing toward high values of ^scaUr a t small neutralino masses. This originates in the 
lower bound on Cg""^ 1 ™ 11 ' reported in Ea. Hll|) . which is effective only for very low neutralino masses (below about 
15 GeV) and becomes more and more stringent as m x decreases. As displayed in Ea. pifl . the size of this lower 
bound, apart from relevant supcrsymmctric details, is determined by the value of (£l x h 2 ) max , which is set here at the 
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value (f2 x /i 2 ) max = 0.3. It is noticeable that at very small values of R, for instance at R — 0.01, all supcrsymmetric 
configurations are within the cosmologically interesting range of Q x (i.e. no configuration of this set is rescaled) and 
provide large values of P^aiax (i.e. large detection rates). 

As we increase the value of R, in our scan we are accessing larger values of m x : again the largest values of C^calar 
are dominated by Higgs-exchange, for Higgs masses close to their lower bound of about 90 GcV. This is also true 
for the annihilation cross section. This approaches its pole at m x ~ therefore, the largest values of ^"^T 011 ^ 

refer to subdominant neutralinos, as m x increases toward m x ~ 45 GeV (which represents the pole in the annihilation 
cross section for the lightest possible A boson). These features are clearly shown in Figs llaljTbl The panel denoted by 
"standard" in Fig llbl refers to the usual case of universal gaugino masses: in this case the neutralino mass is bounded 
from below at about 50 GeV, and therefore all the interesting low neutralino-mass sector is precluded. The last panel 
in Fig llbl (denoted by "global") shows our results for R varied in the interval 0.01 — 0.5: the funnel at low masses and 
the effect of the A-pole in the annihilation cross section are clearly visible. 

We recall that, for each panel at fixed R, the lower value of the neutralino mass is a consequence of the experimental 
bound on the chargino mass, which in turn fixes a lower bound on Mi = R x A/2. The upper value on the neutralino 
mass for each panel is a mere consequence of the fact that we scan the M 2 parameter up to 1 TeV. 

The detailed connection among the values of cr^ 1 ™^ 01 ^ anc ^ tnosc of Q x h 2 is given in Fig. [21 The strong correlation 
between P^aiax ana ^ ^x^ 2 displayed for R = 0.01 reflects the properties of the funnel previously discussed in 
connection with FigEI All the configurations refer to large values of Vl x h 2 : actually, it is the upper bound on the 
neutralino relic abundance which determines the strong bound on the allowed configurations. By changing R from 
0.01 to larger values, we observe that the ensuing increase in m x shifts the configurations of largest ^gcaiax toward 
lower values of relic abundance, as expected from the analitical considerations of the previous Section. From this 
figure we see that a fraction of the largest values of the quantity ^1™!*™' refer to dominant neutralinos, while 
another fraction refers to slightly subdominant neutralinos: 0.01 < £l x h 2 < 0.05. Configurations with fl x h 2 < 0.01, 
even providing the largest values of the scattering cross section (see, for instance, the panel at R = 0.04 in Fig. EJ) 
suffer from a severe rescaling factor £ which somehow reduces their dctcctability. 

The fact that for small values of R the scattering and neutralino-neutralino annihilation cross sections are dominated 
by Higgs exchange is a consequence of two facts: the relatively small values for the lower bounds on mu and ttia 
and the neutralino composition, which, even though dominated by the bino component, nevertheless possesses a non 
negligible higgsino contribution allowing the neutralino to efficiently couple with the Higgs fields. 

Fig. |21 shows that for small values of R (small m x ) the neutralino-neutralino annihilation cross-section is indeed 
dominated by Higgs-exchange diagrams, especially for the largest values of ^"aiar 501 ^ ■ The nrs ^ P an °l of Fig. |3 
which refers to R = 0.01, clearly shows that the annihilation cross section is strongly dominated by Higgs exchange. 
For R = 0.02 the annihilation cross section can be either dominated by Higgs or sfermion exchange: however, the 
configurations which provide values of ^scalar 011 ' m cxccss °f 10~ 8 nbarn (denoted by crosses) show a clear Higgs 
dominance in the annihilation cross section. These features arc progressively lost when R increases: the annihilation 
cross section may be dominated by Z exchange (which, by coincidence, has its pole also at about 45 GeV). 

Finally, Fig. 0] shows that for low values of R, the neutralino composition is dominated by the bino component, 
but a deviation from a pure bino composition is present and is mainly due to a mixture of B with H®. The two 
composition parameters a 2 and a 2 remain aligned along the a 2 + a 2 = 1 diagonal line up to R ~ 0.05, with a clear 
dominance (above 70%) in bino. For larger values of R the correlation between a\ and a 2 starts to deviate from the 
diagonal line, a fact that indicates how the two other components are becoming important (it is mainly 04 which sets 
up). The panel at R = 0.1 shows that the bino component is usually large, but a sizeable mixture starts occurring. 
The last panel in Fig. 0] recalls the situation for the standard case of universal gaugino masses, where the neutralino 
may be any mixture of its component fields. 



IV. CONCLUSIONS 



In the present paper we have focussed our attention on relic neutralinos of light masses: m x < 45 GeV, which are 
allowed in supersymmetric models where no unification of gaugino masses is assumed. We have shown that these 
neutralinos may have elastic cross-sections off nucleons which go up to cr™,^ on ^ ~ 10~ 7 nbarn, with a relic abundance 
of cosmological interest: 0.05 < il x h 2 < 0.3. 
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The present upper limits to £ ^scalar ^ provided by WIMP direct detection experiments 0, IH, do not 

constrain the supcrsymmetric configurations for the light neutralinos considered here. This is especially true once the 
relevant uncertainties (mainly related to the form and parameters of the WIMP galactic distribution function [il| 
and to the quenching factors for bolomctric detectors) arc taken into account. The CDMS upper bound could 
concern a small fraction of supcrsymmetric configurations in the range around 15 GeV, though very marginally, if 
the uncertainties on astrophysical quantities are considered. Moreover, the CDMS bound needs a confirmation by a 
further running in a deep-underground site, as planned by the Collaboration. 

The small-mass neutralino configurations analysed in the present paper are accessible to experiments of direct 
detection with a low-energy threshold and a high sensitivity. An experiment of this type is the DAMA experiment 
with a mass of ~ 100 kg of Nal(Tl), whose results after a 4-years running show an annual-modulation effect at a 4<r 
C.L. which does not appear to be related to any possible source of systematics The DAMA experiment, with its 
high sensitivity, is potentially good to investigate also the relic neutralinos considered in the present paper. 
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FIG. la: Scatter plots of the neutralino-micleon cross section ^"aiar times the rescaling factor £ vs. the neutralino mass, 
for nonuniversal gaugino models with different values of the gaugino mass ratio R = M1/M2: R = 0.01, 0.02, 0.03, 0.04. Crosses 
denote configurations with dominant relic neutralinos (0.05 < Q x h 2 < 0.3), while dots refer to subdominant neutralinos 
(fl x h 2 < 0.05). 




FIG. lb: The same as in Fig. Hal for R = 0.05, 0.1, for the standard value R = 5/3tan 6w — 0.5 and for a generic variation 
of R in the interval 0.01-0.5. 
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FIG. 2: Scatter plots of the neutralino-nucleon cross section (^"alar vs - ^he neutralino relic abundace fi x /i 2 , for R = 
0.01, 0.02, 0.03, 0.04, 0.05 and for the standard value R = 5/3 tan 2 6 W ^ 0.5. 
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FIG. 3: Scatter plots of the fractional amount of the neutralino pair-annihilation cross section due to sfermion exchange vs. 
Higgs exchange, for R = 0.01,0.02,0.03,0.04. Crosses denote configuration for which the neutralino-nucleon scattering cross 
section Scalar * s l ar S er than 10 -8 nbarn. 




FIG. 4: Scatter plots of the neutralino composition in terms of B (ai) and of H° (as) for R = 0.01,0.05,0.1 and for the 
standard value R = 5/3 tan 2 9w — 0.5. The dashed lines denote the line where a\ + a% = 1. 



